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■ Abstract. Travelling waves and conservation laws are studied for a wide class of U{1)- 

CSJ I invariant complex mKdV equations containing the two known integrable generalizations of 

\^ • the ordinary (real) mKdV equation. The main results on travelling waves include deriving 

, ^ I new complex solitary waves and kinks that generalize the well-known mKdV sech and tanh 

solutions. The main results on conservation laws consist of explicitly finding all 1st order 
conserved densities that yield phase-invariant counterparts of the well-known mKdV con- 
served densities for momentum, energy, and Galilean energy, and a new conserved density 
describing the angular twist of complex kink solutions. 



^ ; 1. Introduction 

In this paper, we study travelling wave solutions and conservation laws of general complex 
mKdV-type equations 

ut + auuux + f3u^Ux + ju^xx = (1) 

2 I for u{t,x) with complex coefficients a = ai + ia2, P = /3i + i(32, and real coefficient 7 > 0. 

'sj" ■ An equation of the form ([T]) is equivalent to a coupled nonlinear system 

uit + ((«! + f3i)ul - 2/32M1M2 + (tti - f3i)ul)ui^ 

- ((a2 - f32)ul - 2/3lMlM2 + (^2 + f32)ul)u2x + JUi^xx = 
U2t + (("2 + /32)m? + 2/3iMiM2 + ("2 " f32)ul)ui:„ 

^ . + ((ai - + 2/32M1M2 + (ai + f3i)u2)u2x + lU2xxx = 

■ for the real and imaginary parts of u(t, x) = Ui(t, x) + iu2(t, x). This system reduces to the 

ordinary mKdV equation in the case when U2 = and 02 = P2 = 0. 

Complex mKdV-type equations ([T]) are interesting both physically and mathematically. 
When the coefficients a and (3 are real, such equations describe propagation of short pulses 
in optical fibers [H |2] where the physical meaning of t and x as time and space variables is 
reversed. In the cases where the ratio of these real coefficients is/3/a = 0or/3/a = l/3, the 
resulting equations are integrable systems [31 H], possessing rich mathematical features such 
as soliton solutions which describe nonlinear interactions of two or more travelling waves, 
as well as a hierarchy of conservation laws which involve u, u, Ux, Ux, and increasingly higher 
order x-derivatives of u and u. In contrast, when the coefficients a and /3 are complex, 
little seems to be known about the nature of solutions or the existence of any integrability 
structure for complex mKdV-type equations ([T]). The equivalent coupled systems ©^(E]) 
with ^2 7^ or /32 7^ arise in modelling weakly coupled two-layer fluids 0. 
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Travelling waves for the class of complex mKdV-type equations ([T]) are considered in 
section [2l By use of symmetry reduction and integrating factors, we first derive all smooth 
solutions of the form 

u{t, x) = U{x -ct) = a + hf{x - ct) (4) 
in which a = Oi + m2, 6 = &i + 262 are complex constants, and f{x) is a real valued function 
that either is single-peaked and vanishes for large x (i.e. a solitary wave), or has no peak 
and approaches different constant values for large positive /negative x (i.e. a kink). We next 
derive all smooth solutions having a linear phase 

M(t, x) = exp{i{kx + wt + (f)))f{x — ct) (5) 

where k,w,(f) are real constants, and f{x) is again a real valued function with the same 
general profile (i.e. either a solitary wave or a kink) as considered for the solutions (j4]). Both 
of these two classes of complex travelling waves (jl]) and ([5]) include the familiar sech and 
tanh profiles for f{x) in the case of the ordinary mKdV equation [H]. 

Conservation laws for the class of complex mKdV-type equations ([T]) are then considered 
in section [31 Specifically, by means of multipliers, we derive all conserved densities and fluxes 
of the form 

T(tj X, M, M, Uxj Uxxi Uxx) (6) 
X (t, X, M, M, Uxj Uxi Uxxi Uxxi Uxxxi Uxxxi Uxxxxi Uxxxx) C^) 

which satisfy 

DtT + DxX = (8) 

corresponding to conserved quantities 

/oo 
T dx = const. (9) 
-oo 

for all solutions u{t, x) that have vanishing flux at a; = ±00. The class of conserved densities 
([6]) includes the well-known conservation laws for mass, momentum, energy and Galilean 
energy in the case of the ordinary mKdV equation [7], as well as the first of the higher- 
derivative conservation laws in the hierarchy arising for the two known integrable cases of 
complex mKdV equations. 

In section m the conserved quantities ([9]) are used to explore some features of the travelling 
wave solutions (jl]) and Finally, a summary of the main new results obtained in the 
previous sections is provided in section O 

Hereafter, for convenience we will put 

7 = 1 (10) 
by scaling the time and space variables t — )■ ^/^t, x — )■ y/jx. 

2. Travelling waves 
Each complex mKdV equation in the class ([T]) has the following basic invariance properties: 

scaling x \x,t X^t,u ^ X^^u (11) 

time translation t t + e (12) 

space translation x ^ x + e (13) 

phase rotation u — )■ exp{i(f))u (14) 
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where the parameters are given by A 7^ 0, —00 <e<oo, O<0< 27r. Composition of these 
transformations ( !TT|) -( !Ti|) yields a 4-parameter group of point symmetries admitted by all 
equations ([1]). 

Travelling waves u = U{x — ct) arise naturally as group-invariant solutions [S] with respect 
to the combined space-time translation 

t— j-t + e, X— J-x + ce (15) 

where c is the speed of the wave. From equation ([T]), such solutions satisfy 

- cU' + aUUU' + PU'^U' + U'" = (16) 

which is a 3rd order, nonlinear, complex ODE for U{C,), in terms of the invariant variable 

^ = x-ct. (17) 

Note this ODE f|T6|) inherits invariance under phase rotations 

U exp(i0)f/. (18) 

When the coefficients a and P are real, we can find all real solutions U{^) straightforwardly 
by the use of integrating factors. In contrast, when the coefficients are complex or when we 
seek all complex solutions, the ODE f|T6|) cannot in general be integrated explicitly. However, 
the class of solutions 

U = a + bf{0 (19) 
given by a real function /(^) and complex constants a and b will by-pass these difficulties, 
as we now show. 

2.1. Solitary waves and Kinks. For travelling wave solutions of the form f lT9|) . the ODE 
( !T6l) is given by 

= {adab + (3a% - cb)f + ((a + 2(5)abb + adb^)ff' + (a + f3)b%ff' + bf" (20) 

where, under phase rotations (flSjl . 

a — )■ exp(i0)a, b — > exp(i0)6, f ^ f. (21) 

Consequently, the ODE ( 120|) for /(^) can be simplified by using a phase rotation (12T|) to put 
b = \b\, which then gives 

= {A + Bf + Cf)f + r (22) 

with coefficients 

A = Q;|ap + /3a^ — c, 

B = {a{a + d) + 2l3a)\b\, (23) 
C= (a + /3)|6p. 

Since / is real, these coefficients must be real, and hence we require the conditions 

lmA = lmB = lmC = 0. (24) 

Solving these conditions, we obtain two cases: 

a2 7^ 0, /3i = = /32 = 0; (25) 

as = 0, a2 + (32 = 0. (26) 
For each case, the ODE (l22l) is straightforward to solve by integrating factors. 



We will be interested first in solutions /(^) that describe a solitary wave of the form (jlj) 
having a single peak at ^ = and decaying to a constant value for |^| — )■ oo. Without loss 
of generality, these properties will hold if /(^) satisfies the conditions (i) / — )■ as ^ — )■ ±00 
and (ii) /' = 0, /" 7^ at ^ = ^o- Condition (i) provides asymptotic boundary conditions 

/(±oo) = /'(±oo) = r(±oo) = (27) 

which we impose on the solution of ODE fl221) . Integrating this ODE once, we get 

= Af + \Bf + \Cf + f" (28) 

after the asymptotic boundary conditions ( 1271) are imposed. Then using the integrating 
factor /' and again imposing the asymptotic boundary conditions ( p7|) . we obtain 

= + \Bf + \Cf + f (29) 

which is a separable 1st order ODE. The general solution satisfying condition (ii) is given by 



/(O 



- 6A(:7cosh(v^(e - ^o)) + B' 



exp(±V=A(e-eo)), 
-2B 



± 



A^Q^B = C = Q 
A = 0,fi ^ 

A = B = Q,C 



(30) 



with the coefficients A, B, C given by the two cases fl25|) and fl26l) . This leads to the following 
classification result. 

Proposition 1. Modulo phase rotations a — )• exp(i0)a, b — )■ exp{i(f))b, and translations 
^ — )■ ^ + ^0; travelling wave ODE ^2D^ with asymptotic boundary conditions [2l\ ) has six 
real piecewise- smooth solutions: 



f 



12a/\b\ 
20^2 + 3 



with 



02 + /32 = 0, ai + /3i ^ 



ai + (3i 



, arg 6 = 0; 



(31a) 

(31b) 
(31c) 



/ 



6csch(e)a/|6| 



6 + 4sinh2(e) cosh(sech(e)v^O + 2sinh(e) 



(32a) 



with 



q;2 + ^2 = 0, ai + A 7^ (32b) 



a = J ^tanh(e),arg 6 = (32c) 

Oil + Pi 



where —oo < G < oo; 



with 



where — oo < © < oo; 



with 



where <6 < 27r; 



f^eM-V~c\^\)m 

^2 + /32 = + /3i = 0, ai 7^ or q;2 7^ 
a = sinh(0), arg 6 = 



/ 



12Rea/|6| 



2ccos2(^)^2 + 3 



"2 = A = /92 = 0, tti 7^ 

/ c 

a = \ — exp(i^), arg 6 = 



(33a) 

(33b) 
(33c) 



(34a) 

(34b) 
(34c) 



/ 



6csch(e)|a|/|6| 



^6 + 4cos2(^) sinh2(e) cosh(sccli(e) ^c^) + 2cos(^) siiili(e) 



with 



q;2 = ^1 = /32 = 0, ai 7^ 




tanh(0) exp(i6'), arg 6 = 



where — oo < © < oo and < 9 < 27r; 



(35a) 

(35b) 
(35c) 
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with 



a2 = /5i = /32 = 0,ai^0 (36b) 




, arg 6 = 0. (36c) 

We will next be interested in solutions /(^) that describe a kink of the form (jl]) having 
different asymptotic values for ^ — >■ ±oo and exhibiting an inflection at ^ = ^o- Without loss 
of generality, these properties will hold if /(^) satisfies the conditions (i) / — /± = const, 
with /+ ^ /_ as ^ ^ ±oo and (ii) /" = at ^ = ^o- Since the ODE ([22D for /(O is invariant 
under shifts f f — e, a a + be, b b, we can put /+ = — /_ = /o 7^ by shifting 
f ^ f — (/+ + /_)/2. Hence condition (i) imposes asymptotic boundary conditions 

/(+00) = -/(-oo) = /o, /'(±oo) = r(±oo) = 0. (37) 

Now, integration of ODE by means of the same integrating factors used earlier yields 
the separable 1st order ODE 

= D + Ef + Af + \Bf + \Cf + f (38) 

where D and E are constants. The easiest way we can impose the asymptotic boundary 
conditions fl37|) is by considering the integral curves of this ODE f p8|) in the phase plane 
given by /' = (— -D — Ef — Ap — \Bp — \C f^y/'^ as a function of /, where these conditions 
can be expressed as /'|/=±/o = and \mif^±f^ ^ {1/ f')df = ±00. For these phase plane 
conditions to hold, we must have 



C 

which requires 
and 



^{D + Ef + Af + \Bf + = (/ - UfU + Uf (39) 



B = E = 0,D = lA^/C (40) 



/o = y^Sl/C. (41) 
Combining conditions ( 1iO|) and conditions ( 124|) on the coefficients, we obtain two cases: 

ai = 0, as 7^ 0, f3i = a2 = (52 = 0; (42a) 

ai = 02 = 0, a2 + /32 = 0. (42b) 
In both cases, the general solution of the ODE fl22|) satisfying condition (ii) is given by 



/(O = v/-3A/Ctanh(v^0, A ^ 0, C ^ 0. (43) 



This leads to the following classification result. 
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Proposition 2. Modulo phase rotations and shifts a — )■ exp{i(f)){a + be), b — )■ exp{i(j))b, and 
translations ^ — >■ ^ + ^o? the travelling wave ODE ^2^) with asymptotic boundary conditions 
(37) has two real piecewise- smooth solutions: 

/ = W— sech(e)tanh(sech(e)v/^^0/l&l (44a) 
V a 



with 



"2 = /9i = /32 = 0, ai 7^ (44b) 
^ tanh(9), arg 6 = (44c) 



where — oo < 6 < oo; 



/ = y^^tanh(v^O/l&l (45a) 

with 

02 + /32 = 0, ai + /3i 7^ (45b) 

a = 0, arg 6 = 0. (45c) 

The solutions obtained in Propositions [1] and [2] have interesting analytical behaviour. 
Solutions ([32]), ([35]), ([S]), ([45]) are smooth for all ^, i.e. / e C~ on -oo < ^ < oo. Solution 
fl36]) has a blow-up singularity at ,^ = 0, i.e. |/| — )■ oo as ^ — )■ 0. Solutions fl3T]) and flM]) are 
smooth for all ^ iff c > 0, since otherwise there is a blow-up singularity at some point ^ = 
(such that their denominators vanish). In contrast, solution fl55]) is continuous for all ^ but 
has a cusp at ,^ = 0, i.e. / G C° on — oo < ^ < oo while / G C°° only on < |^| < oo such 
that /' is discontinuous at ^ = 0. 

Theorem 1. For a complex mKdV equation (QP, non-singular solitary wave solutions of the 
form 0) having a single peak at x = at and decaying to a constant value for |a;| — ?■ oo are 
admitted only in the following two cases: 
(i) Ima = 0, (3 = 

Fc I if, , ^ 6sech6 
u{t,x) = e'^i -\ e'^tanhe + 




2 cos 9 sinh 6 + a/6 + 4 cos^ 9 sinh^ cosh.{y/c{x — ct)sech9) 

(46) 

uit, x) = e^\l-{ - , ,, 1 (47) 



a \ 3 + 2c{x — ct)"^ cos^ 9 

where c > 0, a > 0, and the parameters are given by —oo < Q < oo, < 9 < 2tc and 
< 0< 27r. 
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(ii) Im (a + /3) = 

/ X I c / , ^ 6sech0 , 
u{t,x)=e'^J tanhe + , (45 

V « + P V 2 sinh e + V 6 + 4 sinh^ 9 cosh( ^/^(a; - ct)seche) 



<t,x) = e^^/^ (1 - , , , ) (49) 



a + /3 \ 3 + 2c(a; - ct) 

where c> 0, a + P > 0, and the parameters are given by — oo < < oo and < < 27r. 

Theorem 2. For a complex mKdV equation ([7]j, a solitary wave solution of the form ^ 
with a cusp singularity at x = at and a non-singular decaying tail for \x\ oo is admitted 
in the case a + (3 = 0: 

u{t, x) = e''^ (sinh 9 + exp{-^/c\x - ct|)) (50) 
where c > 0, and the parameters are given by — oo < < oo and < < 27r. 
Solutions (jlSD and (jUD generahze the famihar solitary wave 

cosh.{y/c{x — at)) 

which is well known [9j in the case of complex mKdV equations with real coefficients a = a 
and /3 = /3. We note this solitary wave ( 15T!) itself is still a solution in the more general case 
of complex mKdV equations whose coefficients satisfy Im {a + (3) = 0, and it reduces to the 
1-soliton solution in the case of the ordinary mKdV equation (where = 0). 

Solution ( l50l) is somewhat trivial because complex mKdV equations in the case a + /3 = 
reduce to the Airy equation Ut + u^xx = when u has a constant phase. 

Theorem 3. For a complex mKdV equation (CP, non-singular kink solutions of the form 
having an inflection at x = at and approaching different constants for x — ?■ ±oo are admitted 
only in the following two cases: 
(i) Im (a + /3) = 

3c 



M(t,a;) = e^'^W tanh( v/-c/2(x - ct)) (52) 

a + p 



where c < 0, a + /3 < 0, and the parameter is given by < (p < 2-71 . 
(ii) lma = 0, /3 = 



u(t, x) = e 



i<j> 



^ (^i tanh 9 + v^SsechG ta■nh{^y -c/2{x - ct)seche)) (53) 



where c < 0, a < 0, and the parameters are given by — oo < < oo and < < 27r. 



Up to a phase rotation, solution fl52|) is just the familiar kink solution known for the 
ordinary mKdV equation. We see the same kink solution also holds for more general complex 
mKdV equations. Solution f l53|) is a complex generalization of this kink, having asymptotic 
complex values u — )■ e^'^ ^y c/a{itanhQ ± A/3sech6) as x — ±oo. 

Additional features of these solitary wave solutions fH^ - (|3U]) and kink solutions 
will be discussed in section HI 



2.2. Waves with a linear phase. Travelling waves u = U{x — ct) have a natural general- 
ization given by M = exp{i{kx + wt + (f)))U{x — ct) whose form is group-invariant with respect 
to space-time translations combined with phase rotations 

t— )-t + e, X— )-a; + ce,M— )■ exp{i{kc + w)e)u. (54) 

In this generalization the wave is modulated by a linear phase where /c, w are the inverse 
wave length and the frequency of the modulation. A complex mKdV equation ([T]) will admit 
such solutions if f/(0 satisfies the 3rd order, nonlinear, complex ODE 

= i(w - k^)U + i{a - /3)kUU'^ - (c + 3P)U' + aUUU' + /3U'^U' + i3kU" + U'" (55) 

in terms of the invariant variable f|T7j) . This ODE fl55l) cannot in general be integrated 
explicitly. Nevertheless, as we now show, it is possible to find a class of explicit solutions 

U = U = /(O (56) 

where /(^) is a real function. For solutions of this form (156|1 . the complex ODE (!55l) splits 
into a pair of real ODEs 

= f" + ((ai + (3,)f - c - 3e)f' + {(32 - a2)kf (57) 

= 3k f" + (a2 + (32) ff + {w- e)f + (ai - (3i)kf (58) 

constituting an overdetermined system for f{^),k,w. 

This system (I57l) - fl58l) can be solved by an integrability analysis, which we carry out using 
computer algebra, as summarized in section 12.31 The analysis yields all of the different cases 
in which the system reduces to a single ODE for /(^) plus algebraic conditions on k,w. 

Proposition 3. All solutions of the ODE system ^^51^-^5^ for travelling waves with a non- 
zero linear phase are given by: 

f - {c + 3k'^)f + laif^ = (59a) 

with 

a2 = (3i = (32 = 0,ai^0 (59b) 
w = -i3c + 8k^)k; (59c) 

/' T Vc + 3k^f = (60a) 

with 

al + al = (3l + (3^^0 (60b) 

ru = -(3c + 8k^)k (60c) 

± Vc + 3P(ai + /3i) - (as - (32)k = ±Vc + 3k^{a2 + P2) + («i - Pi)k = 0. (60d) 



We will be interested in solutions such that /(^) has the profile of either a solitary wave 
or a kink, as described by the respective asymptotic boundary conditions (1^ or (I57|) . 
Integrating the ODEs (l59l) and ( |60l) under these boundary conditions, we obtain the following 
classification results. 
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Theorem 4. For a complex mKdV equation (|7]j, non-singular solitary wave solutions having 
a non-zero linear phase of the form ^ are admitted only in the case Ima = 0, P = 0: 



uit, x) = \ exp i0 ^ ^ — ' " 61 

V a "^^^ cosh{{x- at) V^+W) ^ ' 

where c > — a > 0, and the parameters are given by —oo < k < oo, < < 271. 

Theorem 5. A solitary wave solution of the form ^ with a cusp singularity at x = at and 
a non-singular decaying tail for \x\ oo is admitted in the case \a\ = 



u{t,x) = Aexp (icp) exp ( -ie^-^^j^(x - _ ^ ct)) exp ( - y ~ ^^D (^2) 

with 0" = Re (a — /3) /Im (a + /3) = — Im (a — f3) /Re (a + /3) 7^ and e = sgn(cr(x — ct)) = ±1, 
where c> 0, > 3, or c < 0, < 3, and the parameters are given by A > 0, < (p < 2tt. 

Solution flM]) is the well-known 1-soliton solution for the integrable case a = a, /3 = in 
the class of complex mKdV equations ([T]). It reduces to the familiar solitary wave solution 
(I5T|) when k = 0. 

In contrast, solution (!62|) is new. By writing it in terms of ^ = x — ct and t, we see that 
its amplitude \u\ = Aexp{—^^/cay{o^'^^3)\^\) has a cusp at ,^ = 0, with an arbitrary height 

A>0, while its phase arg(M) = - e^c/{a^ - 3)^ + e2{cr'^ + l)^c/{a^ - 3ft (modulo 2n) 
exhibits a jump discontinuity at t 7^ yet is continuous at t = 0. Thus, this solution is a 
type of singular peakon. 

Theorem 6. For a complex mKdV equation ([2]j, non-singular kink solutions having a non- 
zero linear phase of the form ^ are admitted only in the case Ima = 0, /3 = 0; 

u{t,x) = + 3fc ) g^p(^^0^ exp{ik{x - (3c + 8P)t)) tanh( v^-(c + 3P) /2{x - ct)) (63) 

V a 

where c < —3k'^, a < 0, and the parameters are given by —00 < k < 00, < < 27r. 

Solution (163|) generalizes the kink solution given by A; = for the integrable case a = 
a, /3 = in the class of complex mKdV equations ([1]). 

Further discussion of the solutions (1611) . fl62|) . (1631) will be given in section HI 

2.3. Computational Remarks. The integrability analysis of the overdetermined ODE sys- 
tem flFn) - fl551) is non-trivial due to the nonlinearity of the two ODEs and the number of 
parameters in each ODE. 

The analysis begins by using differentiation combined with cross-multiplication to elimi- 
nate /"', /", /' until a purely algebraic equation in terms of / is obtained. During elimination 
steps, a case distinction arises whenever the coefficient of a leading derivative can possibly 
vanish. These coefficients involve the parameters ai, a2, Pi, (^2, k, w and sometimes also the 
functions / and /', and so each case distinction increases the size of the system. Because 
the original ODEs (!57j) and (158|) are autonomous, the variable ^ does not appear explicitly 
in the coefficients or the equations. Consequently, any resulting equation that is algebraic 
in / can be split with respect to all powers of /, which gives algebraic conditions just on ai, 
(^2, Pi, P2, k, w. These conditions are first reduced with each other and then used to reduce 
the remaining equations in the system. 
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Because of the constantly increasing number of options of how to proceed in the compu- 
tation (i.e. the number of equations grows due to the repeated case distinctions, offering an 
increasing number of possible reductions to do, which in turn generate new equations), it is 
difficult for a computer program to perform the computations automatically. Computer alge- 
bra programs designed for differential Grobner basis computations (i.e. cross-differentiations 
and cross-reductions of differential equations based on an ordering of derivatives which guar- 
antees the finiteness of the algorithm for linear problems) are typically designed and tuned 
to do a kind of "breadth-first search" that aims to prevent or delay the generation of very 
large equations. But the situation for the present computation is different in that even a 
large algebraic equation involving / is very useful as it can be directly split (because the 
system is autonomous), yielding useful shorter algebraic conditions that involve only ai, 02, 
Pi, P2, k, w. Therefore, it is much more productive to employ a "depth-first search" aiming 
at reducing the differential order of one equation as quickly as possible in the computation. 

For these reasons, the computer algebra package Crack [T3] has been used. The nu- 
merous techniques (modules) in this package for solving overdetermined systems of dif- 
ferential/algebraic equations include automatic and interactive eliminations, substitutions, 
length-shortening of equations, and factorizations, among others. In particular, the ability 
to use the program interactively is essential and gave the complete solution of the integrabil- 
ity analysis. A streamlined version of this interactive run (when executed in an automatic 
mode) used about 300 interactive steps taking 4 seconds in total on a 3GHz PC. In other 
words, the computation is short for a computer program if it is guided properly, but it is 
long enough that a hand computation would be undesirable. (By comparison, the Maple 
program RiffSimp running on a workstation for one day was unable to complete the full 
computation.) 

As a result of tackling this application, the package Crack has been enhanced to perform 
an automatic splitting on equations that contain powers of any function occurring only 
polynomially and having a non-vanishing derivative with respect to at least one variable 
such that this variable neither appears explicitly in the equation nor occurs as a variable of 
another function in the equation. 



3. Conservation Laws 

In the class of complex mKdV equations ([T]), all conservation laws ([H]) can be expressed 
in the integrated form 

lc = -X|+°° (64) 

in terms of corresponding conserved quantities C = T dx, where the conserved density T 
and flux X are functions of t, x, u, u, and x- derivatives of u, u (up to some finite order) after 
eliminating t-derivatives of u, u through equation ([1]). Two conservation laws are equivalent 
if their conserved densities T differ by a total x-derivative -D^T whenever u satisfies equation 
dl]), so that the same conserved quantity C is obtained up to boundary terms T|^°°, where 
T is any function of t, x, u, u, and x-derivatives of u, u. Correspondingly, the fluxes X 
of two equivalent conservation laws differ by a total t-derivative —D{T whenever u satisfies 
equation ([T]). Since any conserved quantity can be decomposed into separately conserved real 
and imaginary parts, there is no loss of generality in considering only real-valued conserved 
densities and fluxes. The set of all such conservation laws up to equivalence for a given 
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complex mKdV equation ([T]) forms a real vector space, on which there is a natural symmetry 
group action generated by scalings (fTTI) . time translations (fT2|) . space translations (fT3|) . and 
phase rotations f lT4|) . 

The ordinary mKdV equation has well-known conservation laws [7J corresponding to con- 
served quantities for mass f^^udx, momentum f^^u^dx, energy /^(3m^ — u'^)dx, and 
Galilean energy f^^^Stu^ — tu^ + xu^) dx. We will be interested, first, in finding all counter- 
parts of these conserved quantities for complex mKdV equations ([!]) through a classification 
of all real-valued conserved densities of the form 

T{t,x,u,u,Ux,Ux) (65) 

modulo total x-derivatives. 

Each conservation law ([8]) of a complex mKdV equation in which the conserved density T 
and flux X are real-valued expressions can be expressed in an equivalent characteristic form 

m 

DtT + Dx{X - r) = 2Re {{ut + auuu^ + (iu^u^ + m...)Q) (66) 
given by a complex multiplier Q = Qi+ iQ2, where the additional flux terms F are linear in 
Ut, Ut and x-derivatives of Ut, Ut (up to a finite order), with u now taken to be an arbitrary 
function oft, x. The characteristic form fl66l) is useful because it leads to a linear determining 
system for all conservation laws (up to equivalence) as follows. 

First, conserved densities T modulo total x-derivatives have a one-to-one correspondence 
to multipliers Q through the variational expressions 

where 5/5u and 5/5u denote variational derivatives (Euler operators) with respect to u and 
u. As a result, Q is a function of t, x, m, m, and x-derivatives of m, u (up to twice the order 
of T modulo D^T). Second, all multipliers Q are determined by a linear system [TOl [TT] 

S S 

—Re {{ut + auuu^ + + m^i:i:)(5) = —Re {{ut + ommm^ + l^u^u^ + Mi.xx)Q) = (68) 

ou ou 

in which u is an arbitrary function of t, x. Since Q has no dependence on ut-, Ut, and 

their x-derivatives, this system ( EE\i has a natural splitting with respect to these variables. 

The resulting split system consists of [121 [IS] the adjoint of the determining equation for 

infinitesimal symmetries, augmented by the Helmholtz integrability equations for variational 

(Euler-Lagrange) expressions. Consequently, the multipliers Q for a given complex mKdV 

equation ([1]) can be characterized as adjoint-symmetries that have a variational form. Third, 

each multiplier Q determines a corresponding conserved density T modulo -D^-T through the 

variational relation flBTj) . In particular, this relation can be explicitly inverted to obtain T 

in terms of Q by means of a homotopy integration formula [TT] 



T = 2Re(M / Q[Xu]dX) + D^T 
Jo 



(69) 



where Q[Am] denotes the function Q evaluated with u, u, and all x-derivatives of u, u replaced 
by Am, Au, and their x-derivatives. 

Now, for classifying conserved densities dHSJ, we will need to find all multipliers of the 
form 

Q(t, X , U, tlx, Uxxj Uxx) ■ ('''0) 
12 



The linear determining system for such multiphers is then given by the determining equation 
for adjoint-symmetries 

DtQ + (2/3 - a)uu^,Q + (a - 2/3)mm^Q + auuD^Q + (3u^D^Q + DIQ = (71) 

and its complex conjugate, augmented by the Helmholtz integrability equations 

B,Re|« =Re|«. Z,J,„^ = I,„|«. (72a) 

UUxx OUx UUxx UUx 

Im = 0, DxRe t: = Re -— , D^.lm -— = 2Im — . (72b) 

OUxx UUxx OUx UUx UU 

In equation f l7T|) . the variables Ut, Ut, and their x-derivatives are eliminated via equation 
(II]) and its x-derivatives (and their complex conjugates). Thus, the system (17T])-(l72l) will 
involve the variables t, x, u, u, and x-derivatives of u, u up to 5th order, with the unknowns 
consisting of Q, a, /3. Since Q is at most 2nd order, this system will split with respect to 
Uxxx^ Uxxx, Uxxxx, Uxxxx, Uxxxxx, Uxxxxx i^to an overdctcrmined system of equations. We use 
computer algebra both to carry out the splitting and to solve the resulting overdetermined 
system, as summarized in section |3TT| which leads to the following classification result. 



Proposition 4. The determining system ((77y-([7^ for 2nd order multipliers ( [7Zy has eight 
solutions: 

Q = 1 (73a) 

with 

ai= 2^1,02 = 2^2-, (73b) 



Q = i 



(74a) 



with 



ai = 2/3i,a2 = 2/32; 



(74b) 



Q = u 



(75a) 



with 



a2 = /32; 



(75b) 



Q 



lU 



(76a) 



with 



ai = 3/3i, a2 = 3/32; 



(76b) 



Q = u 
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(77a) 



with 



with 



with 



with 



ai = 3/3i, a2 = 3/32; (77b) 

Q = iux (78a) 

^2 = ^1 = ^2 = 0; (78b) 

Q = 3uxx + (tti + Pi)u'^u (79a) 

«2 = /32 = 0; (79b) 

Q = StUxx + (oi + Pijtu^u — xu (80a) 

a2 = /32 = 0. (80b) 



The conserved densities flB5]) corresponding to the multiphers fl75]) - fl5n]) can be obtained 
from the integration formula fIBI?]) . One drawback of this formula is that if a multiplier Q 
contains x-derivatives of m or n then the resulting conserved density does not necessarily 
have the lowest possible order. In particular, if Q contains Uxx or Uxx, then T will have 
terms containing Uxx and Uxx, which can be canceled by subtraction of suitable DxT terms. 
Through the Helmholtz integrability equations (!72|) . we can show that these terms are given 
by the integral 

T = 2Re / /' '4M'4 + . r '4M^4) (81) 

J Jo OUxx A io OUxx A 

Hence we obtain an improved homotopy integration formula 

/•I 

T = 2Re / {\uQQ{t,x,\u,\u) — W{t,x,\u,\u,\Ux.i\Ux))— (82) 

JQ A 



where 



W = Ux I {Xux-^;^{t, X, u, u, Xux, Xux, 0, 0) + Xux-^^{t, X, u, u, Xux, Xux, 0, 0) 







ralm ^^(t, X, M, M, Xux, Xux)) dX 



^3) 



in terms of Q\ = Q(t,x,u,u,Ux,Ux, 0,0) and Qo = Qi(t,x,u,u,0,0). Evaluation of this 
formula f l82p for the multipliers fl73l) - fl80l) leads to the following main result. 
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Theorem 7. For a complex mKdV equation (QP, conserved densities of the form / fg3]) are 
admitted only in the following cases: 
(i) Ima = Im/3 

T = uu. (84) 



(ii) Ima = 0,Im/3 = 

T = -?,UxUx + \{a + I3)u^u^; (85) 
T = —StUxUx + h{a + (3)tu^u^ — xuu. (86) 



111 a 



2/3 



T 
T 



u + u; 
i(u — u) 



17) 



(iv) a = 3(3 



T 



(89) 
(90) 



(v) Ima = 0,/3 = 

T = i{uxU — UxU). (91) 

As indicated by their form, the conserved densities (1841) . (1851) and (1861) are phase-invariant 
counterparts of the respective conserved densities for momentum, energy and Galilean energy 
listed earlier for the ordinary (real) mKdV equation. The form of the other phase-invariant 
conserved density fl9T|) suggests that it has the meaning of an angular (phase) twist, since it 
vanishes when the phase of u is constant. The remaining conserved densities flHTl) . fl88l) . fl89|) . 
fl90|) are not phase- invariant. In particular, conserved densities flHTj) and fl88l) are related by a 
phase rotation u — )■ exp(z7r/2)M, while conserved densities fl89|) and fl90l) are similarly related 
by M — exp(z7r/4)M. This suggests that the following complex linear combinations of these 
densities can be viewed as a phase-covariant mass density 

T = u (92) 

and a phase-covariant momentum density 

T = u^ (93) 

where the complex conserved densities (p2l) and (1931) are homogeneous with respect to phase 
rotations (fT4|) on u. Further discussion of the conserved quantities defined by all of the 
phase-invariant densities (15^ . (15^ . (PTj) and phase-covariant densities (1^ . (1551) will 

be given in section HI 

Finally, we will classify all real-valued higher derivative conserved densities having the 
lowest-order form 

T{t, X, It, It, tlx, '^xi '^xxi '^xx) (9^) 

with an essential dependence on Uxx and Uxx, niodulo total x-derivatives. The existence 
of such conserved densities will be important for detecting integrable cases in the class of 
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complex mKdV equations ([T]). For this classification, we need to find all multipliers of the 
form 

Q(ty 11, 11, Uxy Ux: Uxxi Uxxj "^xxxi Uxxxi Uxxxxi ^xxxx) (95) 

with an essential dependence on Uxxxx, Uxxxx, or Uxxx, Uxxx, satisfying the linear determining 
system (1^ . Again, this system splits into the adjoint-symmetry equation (ITTi) and its 
complex conjugate, augmented by the Helmholtz integrability equations which now are given 
by 

2L>^Re = Re (96a) 



XXX 



2Dx\m = Im (96b) 

n OQ dO dO , , 

DlRe - 2DxRe = -2Re (96c) 

l/ f^xxx ^ ^ ^x 

Dll^ - ^Dxlm = -2Im (96d) 
do 

Im ^ = 0, (96e) 



2D,Re = Re (96f) 

^ ^XXXX ^ ^xxx 

SDxlm = 2Im (96g) 



D^Re - 2D,Re ^ = -2Re (96h) 

C/ f^xxx ^ ^XX ^ ^x 

Dllm - 3Dxlm |^ = -6Im (96i) 

OUxx OUx OU 

The determining system ( ITTl) and ( l96l) involves the variables t, x, m, m, and x-derivatives 
of ii, M up to 7th order. Since Q is at most 4th order, this system will split into a large 
overdetermined system of equations, which we derive and solve by using computer algebra, 
as summarized in section ISTTl We thereby obtain the following classification result. 

Proposition 5. The determining system ( [77] ) and for 3rd and 4th order multipliers 
( fg3]) has three solutions: 

Q = iiuxxx + aiuuux) (97a) 

with 

"2 = /3i = /32 = 0; (97b) 



Q = Quxxxx + 2aiix 1x^2^ + ^aiuuuxx + 6ai'Uix^ + 4aiixiXa;'U2: + a^u u (9^ 



"2 = /3i = /32 = 0; 
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(98b) 



Q = ^Uxxxx + QfiiU^Uxx + I4:f}iuuuxx + ^Piuul + l2f}iuUxUx + SPlu^u^ (99a) 

with 

ai = 3/3i, a2 = P2 = 0. (99b) 

To obtain the 2nd order conserved densities (JSj) determined by these muhiphers (|971) - 
(1^^ . we evaluate the integration formula f l69l) and cancel all higher-order terms (containing 
Uxxx, Uxxx, Uxxxx, Uxxxx) by Subtracting suitable total x-derivatives. This yields the following 
classification result. 

Theorem 8. For a complex mKdV equation (QP, conserved densities of the form are 
admitted only in the following cases: 
(i) Ima = 0,/3 = 

T =i\{uxxUx — UxxUx) + iQa{uu^Ux — u^uux); (100) 



2 

T =6uxxUxx - ja{u + u){u- ufiuxx + Uxx) + \a{2uu - - 3m^)(u^ + ul) 



\a{u^ + 14uu + v?)uxUx + la^u'^M^ 



(101) 



(ii) a = 3/3, Im a = Im/3 = 

T =3uxxUxx - f/3(M + m)(m - ufiuxx + M^x) " ^if^i'^'^u^ " 18mm + 27u^){ul + M^.) 
- 1/3(22mm + 3^2 + 3m^)m^m^ + |/3 VmI 



(102) 



The two cases here in which complex mKdV equations admit these higher-derivative con- 
served densities are exactly the well-known integrable cases = and = 3 with real 
coefficients, a = a, /3 = /3. This result indicates that any new integrable cases for com- 
plex mKdV equations having complex coefficients will involve a recursion operator of order 
greater than two. 

3.1. Computational Remarks. For computational purposes, it is convenient to replace 
the complex variables and unknowns by their real and imaginary components, namely, u = 
Ui + iu2, Ux = Uix + iu2x, Q = Qi + iQ2, a = «! + ia2, = A + ^/32- In component form, the 
characteristic form of a conservation law f l66p is given by 

DtT + DxX ={uit + ((«! + f3i)ul - 2(32UiU2 + (ai - l^i)ul)uix - ((«2 - f^2)u\ - 2f3iUiU2 

+ ("2 + P2)ul)u2x + Uixxx)Ql + {U2t + (("2 + ^2^1 + SftWiMa 

+ ("2 - P2)ul)Uix + ((«! - Pl)ul + 2P2U1U2 + («! + Pl)ul)u2x + U2xxx)Q2 

(103) 

in terms of a multiplier pair (Qi,Q2), while the variational expressions (167|) relating such 
multipliers to conserved densities take the form 

Qi = ^, Q2 = ^. (104) 

OUi OU2 

For 2nd order multipliers 

, X, Ml, U2, Uix, U2x, Uixx, U2xx ) (105) 

17 



the determining system (17Ti) - (172l) splits into an overdetermined system of 53 real equations 
for the 6 unknowns <yi,(y2, Pi, P2,Qi,Q2- The main complication in solving the system is 
that it contains terms in which Qi and Q2 appear with ai, a2, Pi, P2 as coefficients, so 
that the system as a whole is nonlinear in the joint unknowns. To formulate and solve this 
system, we use the program Conlaw [U], which in turn calls the package Crack [T3] for 
solving overdetermined systems of differential-algebraic equations. 

Crack uses a wide repertoire of techniques (modules), including eliminations, substi- 
tutions, integrations, direct and indirect separations, length-shortening of equations, and 
factorizations. In the present computation, the typical techniques needed are integration, 
substitution, and separation. Integration of equations is not hindered by the nonlinearity of 
the system, because ai, 02, Pi, P2 are just constants while Qi and Q2 appear only linearly. 
In contrast, indirect separation of equations is sensitive to the nonlinearity. An equation is 
indirectly separable when no unknown in the equation depends on all independent variables 
but each independent variable occurs in at least one unknown in the equation, i.e. direct 
splitting is therefore not possible. Such equations can be separated by a sequence of divisions 
and differentiations, where each division requires a case distinction depending on whether 
the divisor is zero or not. For nonlinear equations, a simple alternative which has been 
implemented in Crack is to differentiate the equation repeatedly and perform case gener- 
ating substitutions and case generating Grobner steps. The entire computation to obtain all 
solutions takes about 1600 steps and is essentially interactive. 

For 4th order multipliers 

, ^1, ^2, ^1x5 '^2x, '^Ixx, '^2xx, '^Ixxx, '^2xxx, ^Ixxxx, '^2xxxx) , 
2(t, X, Ml, U2, Ulx, U2x, Ulxx, U2xx, Uixxx, U2xxx, Uixxxx, U2xxxx) 

the determining system (1711) and (1961) splits into a larger overdetermined system of 212 real 
equations. The resulting system contains, on the one hand, many equations that are linear, 
short (less than 100 terms), and easy to integrate, and on the other hand, equations that 
are longer (up to almost 38000 terms) and nonlinear. Because of this complexity, the whole 
computation to solve this system takes several thousand steps consisting of integrations, 
substitutions, separations, decouplings, and factorizations, while the steps that involve sub- 
stitutions and separations turn out to be very time consuming when performed on the long 
equations. This prompted the following two enhancements to Crack for speeding up the 
computation. 

First, a more specialized separation module has been added as a default module. This 
module performs only direct separations (which are the most frequently needed splittings), 
whereas the general module also performs indirect separations as well as splittings where 
unknowns occur in exponents of independent variables. Splittings with respect to exponents 
are relatively slow, because this involves converting expressions internally between two dif- 
ferent forms (prefix form and standard quotient form). Since none of the unknowns in the 
present computation occur in exponents, a significant speed-up is achieved by avoiding the 
use of the general separation module. 

Second, a new module that partially solves subsystems has been developed, producing an- 
other major speed-up. Previously, each substitution resulting from each integrated equation 
had to be performed in all equations in the system at once, which gets very time consuming 
when the system contains long equations. The new module allows short, linear equations to 
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(106) 



be integrated and/or split, followed by substitution into other short equations, so that after 
several such iterations the equations in the subsystem become much shorter, have low order, 
and are practically all nonlinear in the unknowns. At that stage the preliminary solution is 
then substituted in the whole system including long equations. 

Two other features of Crack greatly aid the present computation. Crack is able to 
collect inequalities (actively generated from conditions on the unknowns in combination 
with algebraic equations that arise during the computation) and use them to select steps 
that will minimize the number of case distinctions and keep the depth of case nestings as 
low as possible. Finally, Crack can be run interactively as well as automatically, which is 
especially useful in solving large nonlinear problems. 

4. Conserved quantities for solitary waves and kinks 

The conserved densities obtained in Theorem [7] for the class of complex mKdV equations 
([1]) yield conservation laws with the following integrated form (IMIl . where 

/CO 
T{t, x,u,u,Ux,Ux) dx (107) 
-oo 

is the integrated density and 

X(t^ X, M, M, Ux, Uxy Uxxi Uxxi Uxxxi "^xxx) (10^) 

is the corresponding flux: 

(i) Ima = Im /3 

/oo 
\u\'^dx = V (momentum) (109a) 

-oo 

X = i(a + /3)mV - UxUx + uuxx + uuxx] (109b) 

(ii) Ima = 0,Im/3 = 

/oo 
-^\ux\^ + \{a + [5)\u\^ dx = 2£ (energy) (110a) 

-oo 

X = "iiuxxUxx - UxUxxx - UxUxxx) + (« + (5){uu'^Uxx + u^uUxx) " (5a + 2(5)uUUxUx 



(110b) 



2 ' 

fOO 



/oo 
— 3t|M^P + |(a + /3)t|u|^ — x\u\'^ dx = 2Q (Galilean energy) (111a) 

-oo 

X = uux + uUx + x{uxUx — uUxx — uUxx) — \x{a + f3)u'^u'^ — t{5a + 2f3)uuUxUx 

- lt{a + 2/3) (m^U^ + ulu^) + 3t{UxxUxx - UxUxxx - UxUxxx) (lUb) 

+ t{a + [5){uxxUu'^ + UxxUU^) + \t{a + (3Yv?v?] 
(iii) a = 2/3 

/oo 
udx = }A (phase-covariant mass) (112a) 

-oo 

X = Pu^u + Uxx\ (112b) 
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(iv) a = 3(3 



/oo _ 
dx = V (phase-covariant momentum) (113a) 

-oo 

X = 2(]u^u-ul + 2uu^^] (113b) 



(v) Ima = 0,(3 = 

/oo 
—i\u\^arg{u)xdx = W (angular twist) (114a) 

-oo 

X = ai{u^uUx — uu^Ux) + iuuxxx — iuuxxx + 2i( ). (114b) 

For a given solution u{t, x) of a complex mKdV equation, these integrated densities fllOQp - 
(11141) will define a finite conserved (i.e. time-independent) quantity (11071) if and only if the 
densities T{t,x,u,u,Ux,Ux) themselves have sufficiently rapid decay as |x| — > oo and also 
the corresponding fluxes (llOSp vanish for x — )■ ±oo. 

We will now determine which of the integrated densities yield a finite conserved quantity 
for the solitary wave solutions (H6l) - (l49l) . (I6T1) and the peakon solution (|62l) obtained in The- 
orems [T] and H] as well as for the kink solutions (1^ obtained in Theorems [3] and O 
To proceed we first write down the asymptotic behaviour of each solution for large ±x such 
that 1^1 > 1/v^- 



Solitary wave ( l46|) : 



Solitary wave (j4 



M = e^('^+^\/-tanhe + 0(exp(-v^|s|)) (115a) 
V a 

lma = 0,(3 = 0. (115b) 



„ = e^(^+^\/- + 0(1/1x1') (116a) 
V a 

lma = 0,(3 = 0. (116b) 



Solitary wave f l48|) : 



u 



-tanhe + 0(exp(-yc|x|)) (117a) 

a + p> 

lm{a + (3)=0. (117b) 



Solitary wave ( I49ll : 



u = e^\^ + 0{l/\x\') (118a) 
\ a + (3 

lm{a + (3) = 0. (118b) 
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Solitary wave ( 16T|) : 



u = 0(exp(-Vc + 3P|a;|)) (119a) 
lma = 0,;g = 0. (119b) 

Peakon ([62]): 



u = 0{exp{-^/cay{a-^ -3)\x\)) (120a) 
|a| = (120b) 

where 

a = Re (a - /3)/Im (a + /3) = -Im (a - (3)/Re {a + (3). (120c) 

Kink (1521): 



M = ie^'^W^^ + 0(exp(-v/^|a;|)) (121a) 
V a + P 

Im(a + /3)=0. (121b) 

Kink (IMD: 



^ = ±e*<^. /^(£±^exp(2A;(x - (3c + 8P)t)) + 0(exp(-v/-2(c + 3F)|x|)) (122a) 
V a 

Ima = 0,/3 = 0. (122b) 
Kink (El: 



M = e^'^±W-(l + 2sech2e) + 0(exp(-v/^|x|)) (123a) 
V a 

Ima = 0,/3 = (123b) 

where 

tan(0± - 0) = ±^ sinh 6. (123c) 

Based on this asymptotic behaviour, it is straightforward to derive the conditions on a 
and P under which each integrated density (I109l) - (lll4p is conserved and finite for each of the 
solutions. The conditions obtained for the solitary wave solutions (HB]) - (H^ . (IMl) and the 
peakon solution (1^^ are listed in Table [T] (where "-" means that an integrated density either 
is defined only for a = /3 = 0, or is infinite, or fails to be conserved). These conditions show 
that all six integrated densities are conserved and finite only for the solitary wave (HHl) in the 
special case = which reduces to the familiar sech solution ( 1511) having the asymptotic 
decay u = 0{exp{—y/c\x\)) for large ±x. For the peakon (162|) . only the momentum (11091) 
is conserved and finite (while both the energy (IllOp and Galilean energy (11111) fail to be 
defined except in the trivial cases a = or l/cr = 0). In contrast, for the kink solutions 
(|52l) and (|53l) . only the angular twist (11140 is conserved and finite. The other kink solution 
(IMj) is exceptional as the linear phase term in its asymptotic behaviour ( 11221) leads to the 
angular twist being neither finite nor conserved. We now list the expressions for all of the 
conserved and finite integrated densities. 
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Proposition 6. Solitary wave l[51\) with speed c > and phase angle < < 27r has the 
following conserved quantities: 



momentum P(c) 



energy £{c) 



a + (3 

6^3/2 



a + (3 

Galilean energy Q = 



phase- covariant mass A4{ 
phase- covariant momentum Vic, (p) 



'a 



(124) 

(125) 
(126) 

(127) 
(128) 

a 

angular twist W = 0. (129) 

Solitary wave ( fglj) with speed c > 0, phase angle < < 2ti, and frequency — oo < — /c(3c + 
8/c^) < oo has the following conserved quantities: 

, 12(c + 3P)i/2 , , 

momentum V{c,k) = , (130) 

a 

6c(c + 3P)V^ 

energy t{c,k) = , (131) 

a 

Galilean energy G = 0, (132) 

angular twist W = 0. (133) 

Interestingly, the momentum V and the energy S for these sohtary waves ( ISTl) and (16T1) 
satisfy the relation 

S = \cV. (134) 
In addition, the Galilean energy is related to the momentum and energy by 

2g = 2t£ - x{t)V = (135) 

where 



1 



oo 



XW = :^ / x\u\'dx (136) 

' J — oo 

defines the center of momentum (in analogy to the definition of center of mass). From 
relations f ll34p and fll35p . this time- dependent quantity ( I136p is given by 

X{t) = at. (137) 

Thus, both of the solitary waves (I5T1) and fl6T|) exhibit the same kinematic relations (1134p 
and fll37p obeyed by a free particle with speed c, whose position x = at coincides with the 
peak of \u\ as determined by ^ = x — at = 0. 

Proposition 7. Solitary waves {"JSj-^J^I ^'^'^ ^^^^ ( f5^) have the conserved angular twist 

W = 0. (138) 

Kink ( f5^) has the conserved angular twist 

W{c, e) = sinh e sech^e. (139) 

a 
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Solitary wave l[61\) has the conserved angular twist 

— '[Oh 

n;(c, k) = Vc2 + 3k^. (140) 

a 

The vanishing of the angular twist for the solitary waves f l46l) - fH^ and the kink f l32]) means 
that the phase angle of u{t,x) is the same (modulo 2-7?) for both x — )■ +00 and x — ?■ —00, 
as shown by the asymptotic behaviour flll5p -( [TT8|) and (1121 p . For the kink ( 153|) . the non- 
zero angular twist (I139P is produced by the finite net rotation in the phase angle of u(t,x) 
between x — )■ +00 and x — )■ —00 as given by the asymptotic behaviour (I123P when 6 7^ 0. 
In contrast, the non-zero angular twist (I140p for the solitary wave (16T]) arises from the linear 
phase k{x — (3c + 8k'^)t). 

Proposition 8. Peakon [6^) with speed c > 0, phase angle < < 27r, and frequency 
—00 < ec'^/2(l — 3cr2)/(cr2 — 3)'^/2 < 00 has the following conserved quantity: 



momentum 




(141) 









Galilean 


phase-covariant 


phase-covariant 


angular 




momentum 


energy 


energy 


mass 


momentum 


twist 


Solitary 
wave (H6l) 


ima = 

/3 = 
6 = 


ima = 

/3 = 
6 = 


ima = 

/3 = 
6 = 






Ima = 

/3 = 


Solitary 
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Table 1. Conditions such that the integrated densities are conserved and finite 



5. Summary 

We have studied travelling wave solutions (|4]) and conservation laws ([6])-(j8]) for a wide 
class of [/(l)-invariant complex mKdV equations ([1]). This class contains the two known 
integrable generalizations of the ordinary (real) mKdV equation given by the Hirota mKdV 
equation and Sasa-Satsuma mKdV equation. Four main results have been obtained. 

Firstly, in Theorems [U to El we respectively derive all solitary waves and all kinks of the 
complex form (jlj), which have the asymptotic behaviour m — )■ a as |x| — )■ 00 and -u — )■ a ± 6/0 
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as X — )■ ±00 where a, b are complex constants and /o is a non-zero real constant. The solitary 
waves in Theorems [1] and [2] describe new complex generalizations of the well-known mKdV 
sech solution, and a new complex rational solution. The kinks in Theorem [3] describe new 
complex generalizations of the familiar mKdV tanh solution. Secondly, in Theorems H] to 
El we respectively derive all solitary waves and kinks having a linear phase of the form (E]), 
with the asymptotic behaviour \u\ — as |x| — t- 00 and \u\ — /o as x — i- ±00. These 
solutions include a new type of complex peakon whose phase has a jump discontinuity where 
the amplitude displays a cusp. 

Thirdly, in Theorem [7| we explicitly find all 1st order conserved densities fl65l) . which yield 
phase-invariant counterparts of the well-known mKdV conserved densities for momentum, 
energy, and Galilean energy, and a new conserved density which describes the angular twist 
of complex kink solutions. Fourthly, we explicitly find all conserved densities fl9^ of lowest 
order in higher x-derivatives. From Theorem [H] the results of this classification establish 
that such higher- derivative conservation laws exist only in the two known integrable cases 
of complex mKdV equations. 

References 

[1] R. Radhakrishnan and M. Lakshmanan, Phys. Rev. E 54 (1996) 2949-2955. 

[2] C. Gilson, J. Hietarinta, J. Nimmo, Y. Ohta, Phys. Rev. E 68 (2003) 016614 (10 pages). 

[3] R. Hirota, /. Math. Phys. 14 (1973) 805-809. 

[4] N. Sasa and J. Satsuma, J. Phys. Soc. Jpn. 60 (1991) 409-417. 

[5] Y. Gao and X.-Y. Tang, Commun. Theor. Phys. 48 (2007) 961-970. 

[6] R. Hirota, J. Phys. Soc. Jpn. 22 (1972) 1456-1458. 

[7] R.M. Miura, C.S. Gardner, M.S. Kruskal, J. Math. Phys. 9 (1968) 1204-1209. 

[8] G. Bluman and S.C. Anco, Symmetry and Integration Methods for Differential Equations, Springer 

Applied Mathematics Series V.154, 2002. 
[9] S.C. Anco, N. Tchegoum Ngatat, M. Willoughby, Physica D 240 (2011) 1378-1394. 
[10] P. Olver, Applications of Lie Groups to Differential Equations, Springer- Verlag, New York, 1986. 
[11] G. Bluman, A. Cheviakov, S.C. Anco, Applications of Symmetry Methods to Partial Differential Equa- 
tions, Springer Applied Mathematics Series V.168, 2010. 
[12] S.C. Anco and G. Bluman, Euro. J. Appl. Math. 13 (2002), 545-566. 
[13] S.C. Anco and G. Bluman, Euro. J. Appl. Math. 13 (2002), 567-585. 

[14] T. Wolf, Applications of CRACK in the Classification of Integrable Systems, in CRM Proceedings and 

Lecture Notes, 37 (2004) 283-300. http://lie.math.brocku.ca/crack/demo/ 
[15] T. Wolf, Crack, LiePDE, ApplySym and ConLaw, section 4.3.5 in: Grabmeier J., Kaltofen E. and 

Weispfenning V. (Eds.), Computer Algebra Handbook, Springer, 2002, 465-468. 



24 



